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ABSTRACT 

We present an analysis of deep HST/WFCS near-IR (NIR) imaging data of the globular 
cluster M4. The best-photometry NIR colour-magnitude diagram (CMD) clearly shows the 
main sequence extending towards the expected end of the Hydrogen-burning limit and going 
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beyond this point towards fainter sources. The white dwarf sequence can be identified. As such, 
this is the deepest NIR CMD of a globular cluster to date. Archival HST optical data were 
used for proper-motion cleaning of the CMD and for distinguishing the white dwarfs (WDs) 
from brown dwarf (BD) candidates. Detection limits in the NIR are around FllQW « 26.5 mag 
and FIQQW « 27 mag, and in the optical around F775IR « 28 mag. Comparing our observed 
CMDs with theoretical models, we conclude that we have reached beyond the H-burning limit 
in our NIR CMD and are probably just above or around this limit in our optical-NIR CMDs. 
Thus, any faint NIR sources that have no optical counterpart are potential BD candidates, since 
the optical data are not deep enough to detect them. We visually inspected the positions of 
NIR sources which are fainter than the H-burning limit in FllQW and for which the optical 
photometry did not return a counterpart. We found in total five sources for which we did not 
get an optical measurement. For four of these five sources, a faint optical counterpart could be 
visually identified, and an upper optical magnitude was estimated. Based on these upper optical 
magnitude limits, we conclude that one source is likely a WD, one source could either be a WD 
or BD candidate, and the remaining two sources agree with being BD candidates. For only one 
source no optical counterpart could be detected, which makes this source a good BD candidate. 
We conclude that we found in total four good BD candidates. 

Subject headings: globular clusters: general — globular clusters: individual(M 4) — stars: brown dwarfs 
— stars: low-mass 


1. Introduction 

Globular clusters (GCs) are the oldest and most 
massive stellar aggregates in our Galaxy. As such, 
they are the best natural laboratories to study 
large, co-eval populations of stars at known dis¬ 
tance and metallicity. Indeed, much of our un¬ 
derstanding of star formation and evolution has 
been derived from observational studies of GGs. 
Nonetheless, we still lack an understanding of the 
very low mass stars (VLMSs) around the faint 
end of the Hydrogen-burning main sequence (MS) 
and of objects beyond that limit, i.e., sub-stellar 
sources, so called “Brown Dwarfs” (BDs). This is 
especially true for objects at the low-metallicities 
typical for the GGs in our Milky Way (see below). 

BDs present a link between stars and plan¬ 
ets, and thus are important for our understanding 
of both star and planet formation and evolution. 
BDs are sub-stellar objects that are not massive 
enough to ignite and sustain Hydrogen burning. 
Thus, where low-mass stars will retain their lumi¬ 
nosity for a Hubble time or longer, a BD will con¬ 
tinue to cool and become fainter with age. Like 
giant gas planets, BDs have complex atmospheres 
([Burrows et al. 19W1 IBurrows et al. 200T1l . The 
distinction between stars, BDs and planets is ei¬ 
ther based on mass or on formation. In gen¬ 
eral, stars have masses > 80Mj and can sustain 
Hydrogen burning, BDs have masses between 80 


and 13 Mj and cannot sustain Hydrogen burning 
but a short period of Deuterium burning, and gi¬ 
ant planets have masses below 13 Mj and cannot 
sustain Deuterium burning (Burrows et al. 1997, 
Stamatellos 2014, but see also Sect. 3.3). Be¬ 


cause low-mass stars and BDs have life-times much 
longer than the age of the Galaxy, the GG VLMSs 
and BDs are also important tracers of Galactic 
formation and chemical evolution. 

The formation of BDs is a matter of con¬ 
siderable dispute. They might have formed in 
the same way as (low-mass) stars from tur¬ 
bulent cloud fragmentation (Elmegreen 1999 
IWhitworth fc Goodwin 200^ lAndre et al. 2012lb 
which would imply a continuous extension of the 
IMF into the sub-stellar regime. On the other 
hand, BDs might form e.g. from the ejection of 
stellar embryos or sub-stellar clumps which did 
not have the chance to accumulate enough mass 


( Reipurth fc Clarke 2001[[Basu &: Vorobyov 2012 ) 
Kroupa & Bouvier (2003) suggested that BDs 
form via photo-evaporation of protostars through 
nearby massive stars. This might also suggest 
an increase in the number of BDs with cluster 
mass, as more massive clusters have more O stars 
which can produce more BDs. BDs might also 
form from the fragmentation of circumstellar disks 
(Stamatellos et al. 2011, Thies et al. 2010, Kaplan 
et al. 2012). In this case, the number of BDs in 


2 















clusters could be enhanced in dense clusters as 
dynamical interactions between cluster stars lead 
to more disk fragmentation. Thies et al. (2015) 
concluded that BDs likely form not just via one 
formation scenario, but from a combination of 
various channels. 

Large surveys undertaken in the past decade 
have detected large numbers of BDs. For example, 
the Two Micron All Sky Survey (2MASS; Skrut- 
skie et al. 2006), the Sloan Digital Sky Survey 
(SDSS; York et al. 2000), the United Kingdom 
Infrared telescope Deep Sky Survey (UKIDSS; 
Lawrence et al. 2007), the Wide-held Infrared 
Survey Explorer (WISE, Wright et al. 2010) all 
sky survey have been very successful in hnd- 
ing such cool, low-mass objects]^ However, for 
most of these BDs, key physical properties like 
metallicity and age are unconstrained (see e.g. 
DwarfArchives.org). In fact, determining the 
physical parameters is extremely difficult and a 
major hurdle in BD research. Observations of 
open star clusters and star forming regions, where 
all sources are at the same distance and metal¬ 
licity, can mitigate this problem (e.g. Steele et 
al. 1995, Rebolo et al. 1996, Martin et al. 2001, 
Pinheld et al. 2003, Boudreault & Lodieu 2013, 
Casewell et al. 2014), but only for young and 
metal-rich objects. Thus the need for benchmark 
sources is especially evident for the metal-poor 
regime, and indeed we still do not know much 
about old, metal-poor BDs. So far, only very few 
low-metallicity, old VLMSs near the H-burning 
limit, and even fewer sub-stellar (candidate) halo 
objects have been identihed (e.g. Burgasser et 
al. 2003, 2009, Lepine et al. 2004, Burgasser & 
Kirkpatrick 2006, Cushing et al. 2009, Sivarani 
et al. 2009, Murray et al. 2011, Mace et al. 2013, 
Pinheld et al. 2014, Luhman & Sheppard 2014, 
Burningham et al. 2014, Kirkpatrick et al. 2014). 

This is where GCs come in: they are massive, 
and thus might have produced BDs in large num¬ 
bers, and they are also the oldest and most metal- 
poor stellar aggregates in our Galaxy. Potentially, 


^Compilations of BDs and low-mass stars can 
be found on DwarfArchives.org, the “List of 
Brown Dwarfs” maintained by W. R. Johnston on 
http: //www .johnstonsarchive.net/astro/browndwarflist.html 
and the “List of all ultracool dwarfs” maintained by J. 
Gagne on https://jgagneastro.wordpress.com/list-of- 
ultracool-dwarfs / 


GCs are the ideal hunting ground for old, metal- 
poor benchmark VLMSs and BDs which are much 
needed if we are to test stellar and sub-stellar 
formation and evolution theories and models of 
metal-poor (sub-)stellar atmospheres. 

However, identifying substellar objects in GC is 
challenging due to their intrinsic faintness. There¬ 
fore, the closest GCs make the best targets for this 
kind of research. Out of the GCs in our Galaxy, 
M4 (NGC6121) is the closest GG to us; distance 
estimates range from 1.7 kpc (Hansen et al. 2004) 
to « 2 kpc (e.g. Bedin et al. 2009). Braga et al. 
(2015) estimated a true distance modulus of 11.28 
mag based on RR Lyrae period-luminosity and 
period-Wesenheit relations, resulting in a distance 
of 1.8 kpc, which agrees well with previous esti¬ 
mates. Malavolta et al. (2014) analyzed 7250 spec¬ 
tra for 2771 cluster stars and found a metallicity 
of [Fe/H] = -1.07 dex (RGB) to [Fe/H] = -1.16 
dex (sub giant branch and MS stars), which is well 
in line with previous metallicity estimates (e.g. 
[Fe/H] = —1.15 dex according to the 2010 update 
of the Harris, 1996, catalogue of globular clusters 
in the Milky Way). 

As such, M4 is a prime target for ultra-deep 
observational studies. Indeed, deep optical stud¬ 
ies with the Hubble Space Telescope (HST) have 
been undertaken by Richer et al. (1997, 2004) and 
Bedin et al. (2009), yielding impressive results. 
Richer et al. (1997, 2004) estimated the fraction 
of similar-mass photometric binaries to be small 
(just 2% in their outer field, falling to just 1% to¬ 
wards the cluster core). Milone et al. (2012) sug¬ 
gested a much higher total binary fraction, raising 
from 10% at the halfmass radius to 15% towards 
the cluster core. The present day mass function of 
the lower-mass MS stars is flat (Bedin et al. 2001), 
with a slope of a = 0.1 and a further flattening to¬ 
wards the cluster centre (Richer et al. 2004). The 
cluster WDs suggest that the initial mass func¬ 
tion (IMF) above 0.8 Mq was much steeper than 
the present day mass function. Bedin et al. (2009) 
presented the deepest optical colour-magnitude di¬ 
agram (GMD) to date of this cluster and located 
the faint end of the WD cooling sequence in M 4 at 
F6061U = 28.5 mag, suggesting an age of 11.6±0.6 
Gyr. This agrees with the finding of Hansen et 
al. (2004), who found a WD based age of 12.1 
Gyr. The ongoing HST M4 core project (PI L. 
Bedin, GO-12911) searches for binary dark com- 
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panions to MS stars. The high-accuracy astrome¬ 
try and photometry of this data-set, together with 
archival material and a part of the deep near-IR 
(NIR) data presented in this paper, have already 
been used to identify two distinct sequences along 
the lower-mass MS (|Milone et al. 2014)1 . 

M4 is also so far the only GC known to host 
a planetary system, PSR B1620-26 (Sigurdsson et 
al. 2003), which challenged the planet-metallicity 
relation of the standard planet formation model 
at that time (|Fischer fc Valenti 2005)1 . Recently, 
Hasegawa & Hirashita (2014) suggested that the 
critical metallicity for gas giant formation is 
[Fc/HJr:! 1.2 dex, which agrees with M4’s metal¬ 
licity. Beer et al. (2004) suggested a metallicity- 
independent formation scenario, in which the 
planet in M 4 formed through dynamically induced 
instability in a circumbinary disc. If true, then we 
can expect many planets to form especially in 
the dense GCs in which dynamical interactions 
between cluster stars are ubiquitous. Since BDs 
might form in a similar way, we then might also 
expect that many more BDs form in dense GCs 
compared to open clusters. 

This paper is structured as follows. In Sec¬ 
tion [2] we describe the observations and reduction 
of our NIR and archival optical data. In Section 
we present and discuss the NIR and optical-NIR 
CMDs, and our results and conclusions are sum¬ 
marized in Section m 

2. Observations 
2.1. NIR Data 

The NIR observations of the globular cluster 
M4 were carried out in April 2012 with the Wide 
Field Camera 3 (WFC3) on board the HST, us¬ 
ing the FWQW and F160IF filters (program GO- 
12602, PI: Dieball). All observations were made 
at a single pointing position on a field centred 
at about one core radius North East from the 
cluster centre. This region in the cluster had 
been the focus of the programs GO-5461 (Richer 
et al. 1997, 2004) and GO-10146 (Bedin et al. 
2009) and thus is fully covered by deep opti¬ 
cal observations. A standard 4-point WFC3- 
IRDITHER-BOX-MIN dither pattern and a sam¬ 
pling of NS AMP 14 SPARS50 was applied dur¬ 
ing the observations to get a well-sampled point 
spread function (PSF). The WPC3 field of view 


is 136" X 123" in the IR channel, with a resolu¬ 
tion of 0.13" X 0.121" per pixel. All NIR obser¬ 
vations were carried out during two consecutive 
orbits on April 16th (FllOW) and four orbits on 
April 20th to 21st (F1601F), comprising a total of 
8 individual images in FllQW and 16 images in 
F1601F, each 653 seconds, resulting in total ex¬ 
posure times of 5223 seconds (FllOlF) and 10447 
seconds (F1601F). 

Using the pipeline produced flat-fielded (FLT) 
images, we first created a master image for each 
of our IR filters, using multidrizzle running un¬ 
der PyRAF (the Python-based interface to iraiQ. 
multidrizzle corrects geometric distortions that 
are present in the input images and combines them 
into a master image. Shifts between the individ¬ 
ual images are expected, as we have applied a 
dither pattern. On top of this, telescope breath¬ 
ing can affect guide star tracking and as a result 
can cause small shifts, typically on a sub-pixel 
scale (see e.g. the multidrizzle handbook avail¬ 
able on the STScI webpages). In order to ascertain 
that all shifts are taken into account, we created 
geometric-distortion corrected individual images. 
Based on the coordinates of the same 10 stars 
in each of these images (selected to be well dis¬ 
tributed over the field of view), accurate shifts be¬ 
tween the individual and the master images were 
then determined using tweakshift. 


The master images created in this way are dis¬ 
played in Figures and Note that these master 
images serve as reference images for the positions 
of the stars detected by DOLPHDT, but the pho¬ 
tometry was actually performed on the individual 
images (see 2.2 below). 


2.2. Photometry of the NIR data 

Photometry was performed on the individual 
FLT images using the DDLPHQT software]^ devel¬ 
oped by A. Dolphin as a generalized version of 
HSTphot (Dolphin 20001. DDLPHQT runs on the 
FLT images downloaded from the STScI archive, 
i.e., no further processing of the images is neces¬ 
sary nor recommended. Photometry is performed 


^IRAF (Image Reduction and Analysis Facility) is dis¬ 
tributed by the National Astronomy and Optical Observa¬ 
tory, which is operated by AURA, Inc., under cooperative 
agreement with the National Science Foundation. 

®http: //americano.dolphinsim.com/dolphot/ 
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on each of the input images, using the WFC3 mod¬ 
ule which replaces the analytic PSF model with 
a look-up table computed using Tiny Tim PSFs 
(Krist & Hook 2011). For our data, we used Jay 
Anderson’s PSF libraries for the WFC3/IR fil¬ 
ters and the pixel area maps available from the 
DOLPHOT webpage. The WFC3 module also in¬ 
cludes a photometric calibration to the VEGA- 
mag system. As a reference frame for a common 
physical (image) coordinate system we use our 
deepest image, i.e., the master image created with 
multidrizzle. The DOLPHOT package also pro¬ 
vides routines to mask all pixels that are flagged 
as bad in the data quality arrays, to multiply with 
the pixel area maps, to calculate sky images, and 
to align all input images (or rather the source coor¬ 
dinates found in each input image) to the reference 
frame (our drizzled master images) using user de¬ 
fined lists of stellar coordinates in each input and 
the reference image. For these steps, we used the 
recommended WFC3 IR parameter settings. 

Note that performing photometry on the driz¬ 
zled image is expected to provide sub-optimal pho¬ 
tometry, because the drizzling process affects the 
PSF and the noise characteristics of the drizzled 
images. Instead, DOLPHOT runs on all input im¬ 
ages simultaneously and thus is capable of provid¬ 
ing deep photometry. We started with the recom¬ 
mended parameters, and then refined some param¬ 
eters to push our photometry as deep as possible]^ 

DOLPHOT can run on data from multiple filters 
and cameras, thus we were able to do the pho¬ 
tometry on both FllOkF and F160IF data sets 
simultaneously. The output file includes x and y 
positions, photometric parameters like sharpness, 
crowding (the magnitude difference to the mea¬ 
sured magnitude if no neighbouring sources would 
be fit simultaneously), the object type, and mag¬ 
nitudes in the VEGAmag system for all sources 
detected in the FllQW and FIQQW data. A cal¬ 
ibrated NIR CMD is plotted in Fig. As can 
be seen, the GMD is exceptionally deep. The to¬ 
tal catalogue contained 51 311 sources, but a large 
number of those will be spurious detections. 


■^We used Forcel = 1, FlagMask = 4 to eliminate satu¬ 
rated stars, WFCSIRpsfType = 1 for the Anderson PSF 
cores, FitSky = 2, SigFind = 1.5, SlgFindMult = 0.8, 
SigFinal = 1.5, and RPSF = 15. 


2.3. Optical Data and photometry 

Among the archival HST material the only im¬ 
ages which offer any hope of detecting BDs in the 
optical are those deep and in the reddest filters 
available, i.e., FllhW and F8141E (indeed, tech¬ 
nically >2/3 of these pass-bands are already in the 
NIR wavelength range). 

We extensively searched in the FIST archive, 
and identified four ~1200 seconds deep images 
taken with ACS/WFG in FllbW under program 
GO-10146 (PI: Bedin) as the optimal for our pur¬ 
poses, as those images were taken in low-sky mode, 
are well dithered, and collected in a well defined 
epoch in 2005.48 (i.e., about 6.8 years before the 
GO-12602 data). 

Images were reduced with the software de¬ 
scribed in great detail in Anderson et al. (2008). 
Briefly, the method is essentially a PSF-fitting 
where all pixels from all images are simultaneously 
fitted using the appropriate PSFs which account 
for the spatial and temporal shape-variation on 
each individual image. The key to the method 
is an optimal knowledge on how to transform 
those pixels into a common reference frame and an 
exquisite empirical PSF modeling. This software 
is well tested and used in all the twelve works of 
the series “The AGS Survey of Galactic Globular 
Glusters” (see, Sarajedini et al. 2007). 

The first photometric run returned 19990 opti¬ 
cal detections. As the detection of BDs might be 
only marginal in these optical images, we relaxed 
the finding criteria in a second photometric run, 
imposing to save all the local maxima, and not 
only the significant ones, resulting in 1.5 million 
detections. Although most of these local maxima 
likely are just fluctuations of the background noise, 
we can still use them to set an optical upper limit 
to the NIR-detected BD candidates. 

For visual inspection of the WDs and BD candi¬ 
dates, stacked images were created from the CTE- 
corrected FLT images. The stacked images have 
pixels values resulting from median clipping of the 
GTE-corrected FLT individual images, and are su- 
persampled by a factor of 2 in each direction (see 
Anderson et al. 2008 and Anderson & Bedin 2010 
for further details). 
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3. The Colour-Magnitude Diagrams 


We used the data in the two NIR filters FllQW 
and F160W to create a NIR CMD. As the full NIR 
catalogue contains a large number of spurious de¬ 
tections, we created a best-photometry NIR CMD 
(described below in Sect. 3.1). However, since 
WDs and BD candidates have similar magnitudes 
and colours in the NIR, BD candidates cannot be 
distinguished based on the NIR data alone. We 
then used additional optical F775W data to cre¬ 
ate proper motion cleaned optical-NIR CMDs (see 
Sect. 3.1) in which the WD and MS sequence are 
clearly separate. Thus, the optical data are used 
to distinguish the WDs from the BD candidates. 


3.1. Best-photometry NIR CMDs 


In order to produce a clean CMD of only “good” 
stellar sources, we selected the output catalogue in 
such a way that the faint sources will be dominated 
by true stellar detections while keeping the num¬ 
ber of spurious detections at bayj^ The resulting 
best-photometry catalogue includes 2526 sources, 
i.e, « 5% of all detections. The best-photometry 
CMD is remarkably clean, as can be seen in Fig.|^ 
We can clearly see the MS delineating down to¬ 
wards the expected H-burning limit and going be¬ 
yond that to fainter sources. Note that this CMD 
is not proper motion cleaned, as we show all NIR 
sources that satisfy our selection criterion on the 
photometry. All sources with magnitudes fainter 
than 24 mag in FllOW were visually inspected on 
the FllOW master image, resulting in 177 visually 
confirmed NIR sources around or fainter than the 
H-burning limit in FllQW (see Sect 3.3 below). 


The MS is narrow between the two MS “knees” 
{Film « 15 mag and FlIOW « 18 mag). The 
first “knee” occurs around T^ff Ki 4500 K and a 
mass of « O.55M0 and is due to the formation 


®For this purpose, we settled on a selection that only in¬ 
cludes source with object type 2 or less (i.e., only “stars”). 
We selected a sharpness between -0.05 and 0.05 (a sharp¬ 
ness of zero denotes a perfectly-fit star, a negative value 
indicates a broader source, a positive value indicates a 
source too “sharp”, for example a cosmic ray - for an un¬ 
crowded field, sharpness values between -0.3 and 0.3 are 
recommended in the DQLPHOT manual, so we apply a stricter 
selection criterion here). The crowding parameter indicates 
how much brighter a source would be if nearby stars would 
not have been measured simultaneously. A crowding of zero 
indicates an isolated star. We allowed for a crowding of no 
more than 0.3. 


of molecules in the cool stellar atmospheres, and 
the second knee is a result of increasing electron 
degeneracy in the (sub)-stellar interior close to the 
H-burning limit (Baraffe et al. 1997). Below the 
second knee, the MS broadens towards fainter and 
lower-mass MS stars. Indeed, the low-mass MS 
splits into two branches, as shown in Milone et al. 
(2014), who had first pointed out multiple stellar 
generations among VLMSs in M 4, based on high- 
precision deep optical HST data from the HST 
M4 core project and our NIR data. The split can 
be clearly seen in the NIR CMD (see Milone et 
al. 2014, their Fig. 2), as opposed to the optical 
data, demonstrating that the NIR is also an ideal 
waveband to search for multiple sequences along 
the lowest-mass and hence faintest MS. Note that 
the goal in Milone et al. (2014) was to search for 
multiple generations along the low mass MS based 
on very precise photometry. In contrast, in this 
paper, our goal is to go deep and well beyond the 
H-burning limit. 

In our Fig. we can see the MS delineating 
further down towards the expected end of the H- 
burning sequence, marked with red slashed lines 
and light-red shaded area. See Sect. |3.3| below for 
a discussion on the mass and NIR magnitude at 
the H-burning limit. Around FllQW « 24.5 mag 
the number of sources on the MS decreases and 
the MS peters out. On the blue side of the faint 
MS, the WD sequence can be seen, starting around 
FllOW < 22 mag and FllOW - F160W « 0.2 
mag and going fainter and redder. 


An increase in source number seems to be ap¬ 
parent below FllOW > 25 mag, i.e., below the 
expected end of the H-burning sequence. This is 
the area in the CMD where we expect the BDs to 
appear. This happens to coincide with the WD 
sequence as well, i.e., this is the region where WD 
and BD cooling sequences would be expected to 
cross. Unfortunately, this also means that we can¬ 
not disentangle WDs and BD candidates based 
on our NIR data alone. In order to help with 
both cluster membership determination and dis¬ 
tinguishing WDs and BD candidates, the deep op¬ 


tical data from GO-10146 were used (see Sect. 2.3 
and Sect 


3.2). 


The WD as well as the BD regions have been 
indicated in Fig. For orientation purposes, we 
have plotted a WD sequence (blue line). The WD 
cooling sequence was constructed by interpolating 
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on the Wood (1995) grid of theoretical WD cooling 
curves, adopting a mean WD mass of 0.55Mq. Us¬ 
ing a grid of synthetic DA WD spectra kindly pro¬ 
vided by B. Gansicke (see Gansicke et al. 1995) we 
carried out synthetic photometry with PySynphot. 
Note that we have shifted the WD cooling se¬ 
quence to get a reasonable match to the under¬ 
lying GMD. This required a rather large distance 
of 2.2 kpc and a reddening of E{B—V) = 0.55 mag 
(a standard reddening law (Seaton 1979) is built 
into PySynphot). Our cooling sequence starts at 
Teff = 50 000 K and terminates at Tg// = 8 000, 
but note that the coolest WDs in M 4 have temper¬ 
atures as low as T^ff = 4000 (Bedin et al. 2009). 

In addition, we have marked the location of 
the known field BD SDSS-J125637.13-022452.4 
dBurgasser et al. 2009D , which has a metallicity 
similar to M 4 but is likely several Gyr younger. As 
a consequence, its cooling time is shorter and thus 
it is expected to be brighter than the M4 BDs. 
The observed J- and i7-band magnitudes agree 
with a 5 Gyr old, O.O78M0 source at a metallic¬ 
ity of [M/H]=-0.5 dex. We scaled the observed J- 
and i7-band magnitudes to the WFG3 NIR filters 
and applied M 4’s distance and reddening (follow¬ 
ing Hendricks et al. 2012, we adopted a reddening 
oi E{B — V) = 0.37 mag, a true distance modu¬ 
lus of 11.28 mag, R = 3.67, Aj = 0.302 * Ay and 
Ah — 0.191 * Ay) and over-plotted the field BD 
on our GMD. Its location supports that our data 
are indeed deep enough to reach well into the BD 
zone. 

Two 12 Gyr isochrones have been overplotted 
on the best-photometry GMD, a BT-Settl model 
based on the Asplund et al. (2009) solar abun¬ 
dances and the Barber & Tennyson (2006) line 
list (Allard et al. 2012), for a metallicity of 
[M/H] = —1 dex; and a Dartmouth model for 
[Fe/H] = —1.2 and [a/Fe] = -1-0.4 (Dotter et al. 
2008). Note that we do not attempt to derive clus¬ 
ter parameters from the isochrone fitting, instead, 
we have fit the isochrones by eye so that they best 
overlap with the underlying GMEj^ Note also that 

®The best fit was achieved with the reddening law from Hen- 
dricks et al. 2012 but a larger Ay = 1.9 for the Dart¬ 
mouth isochrone. For the BT-Settl isochrone, we chose a 
smaller distance modulus of 11 mag and Ay = 1.8. The 
difference in shape of the isochrones, as well as the differ¬ 
ence in the best-fit parameters, reflect the differences in 
the underlying physics, i.e. treatment of the stellar atmo- 


both sets terminate at a stellar mass of 0.083 Mq 
(BT-Settl) or 0.1 Mq (Dartmouth), i.e., they do 
not reach to the H-burning limit. 


3.2. Optical-NIR CMDs 


Our NIR and optical catalogues have been 
matched using a six-parameter linear transfor¬ 
mations between the star positions in the dif¬ 
ferent epochs. As reference stars for the trans¬ 
formations we used only well-measured, isolated, 
non-saturated cluster stars with a high signal-to- 
noise and low residuals. The predicted positions 
of the first epoch sources in the second epoch are 
compared with the observed positions and the dis¬ 
placements between first and second epoch are 
calculated. The top panel in Figures and 

12 shows the displacements in WFC3/IR, based 
on the total DOLPHOT NIR (51 311 detections) and 
optical (19990 sources) catalogue. Since cluster 
stars have been used for the reference list, we ex¬ 
pect cluster members to agglomerate around zero 
in the displacement vector point diagram. Indeed, 
two populations can be distinguished: a dense and 
tight agglomeration of data points around AY = 0 
and AX = 0 which denotes the cluster members, 
and a more widely spread data region centering 
around AY « 0.75 and AX « —0.5 which denotes 
field sources. The latter are mostly Bulge sources, 
reflecting the low tangential motion of M4 around 
the Galactic center (see also Bedin et al. 2003). 

The corresponding GMDs are plotted in the 
second row in Figures Hill] and The left 
GMDs show all sources within up to two WFG3 
pixels displacement, the middle CMDs show only 
sources with a displacement of no more than 0.1 
pixels which suggests that they are cluster mem¬ 
bers, and all non-cluster sources are shown in the 
right diagrams. A displacement of 0.1 pixels corre¬ 
sponds to a proper motion of 1.9 mas/year, based 
on our timeline of 6.8 years and a pixel scale of 
0.13". This agrees well with e.g. Bedin et al. (2003, 
2009) and Zloczewski et al. (2012). As can be seen, 
the proper-motion cleaned “cluster” CMDs (mid¬ 
dle diagrams) show a well defined MS, terminat¬ 
ing around F775W « 26 mag in both CMDs in 


spheres including molecules. For a more in-depth discus¬ 
sion of the input physics to the models we refer the reader 
to the BT-Settl and Dartmouth webpages and the refer¬ 
ences given there (https://phoenix.ens-lyon.fr/Grids/BT- 
Settl/ and http://stellar.dartmouth.edu/models/). 
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Figures [5] and as well as a well defined WD se¬ 
quence (light-blue data points) going down to the 
bottom of the WD sequence around FllbW « 28 
mag. 

On the other hand, the proper-motion cleaned 
optical-NIR CMDs do not show any BD candi¬ 
dates, i.e., sources fainter than FllbW > 28 
mag and redder than F775W — FllQW > 3.5 or 
F775W-F160W > 4.5 mag (see also Fig.0. This 
was expected, as the deep optical CMD presented 
in Bedin et al. (2009) did not show any potential 
MS sources fainter than the H-burning limit. How¬ 
ever, and most importantly in the context of this 
paper, the BDs are expected to be much redder 
than the WDs in the optical-NIR CMDs. Indeed, 
the WD sequence and the MS are clearly separated 
in our optical-NIR CMDs. Therefore we can use 
the deep optical data set to identify the WDs in 
our NIR CMDs and disentangle WDs from poten¬ 
tial BDs. 

3.3. The minimum mass at the Hydrogen¬ 
burning limit 

What is the minimum mass at the Hydrogen¬ 
burning limit in a low-metallicity cluster like M 4? 
And, as a consequence, at which NIR magnitude 
and colour do we expect the H-burning limit in 
our NIR CMD? Early theoretical work (Kumar 
1963) suggested a lower limit to the stellar MS 
at « O.O 7 M 0 for population I and « O.O 9 M 0 for 
the metal-poor population H stars, to which also 
GCs belong. Hayashi & Nakano (1963) suggested 
that stars less massive than O.O 8 M 0 cannot un¬ 
dergo Hydrogen burning. They further suggested 
that the limiting mass is not very different for 
Population I and Population H stars. Burrows et 
al. (1993) presented a zero metallicity theoretical 
model (i.e.. Population HI) which suggests a lim¬ 
iting H-burning mass as high as O.O 94 M 0 . Treat¬ 
ment of the atmosphere has a considerable impact 
on the predicted limiting mass, as can be seen in 
Chabrier et al. (2000) who suggested a limiting 
mass of O.O65M0 for models that include dust for¬ 
mation, Saumon & Marley (2008) who suggested 
a limiting mass of O.O 75 M 0 for cloudless mod¬ 
els, and O.O 7 OM 0 for cloudy models; all for solar 
metallicity and an age of 10 Gyr. 

Previous theoretical works suggest that the H- 
burning limit is at higher masses for more metal- 
poor stars. The models roughly agree on a limiting 


mass of O.O 75 M 0 for solar metallicities. Following 
Hayashi & Nakano (1963), we conservatively as¬ 
sume that the H-burning mass for a population 
as metal-poor as M4 is between O.O 75 M 0 and 
O.O 8 M 0 . 

Unfortunately, detailed sub-stellar models for 
sub-solar metallicities are presently not available. 
Updated models that extend well into the BD 
regime are currently being computed (Allard, pri¬ 
vate communication). 

However, the most recent set of BT-Settl mod¬ 
els (Allard et al. 2012 ; 2013) suggest a strong 
metallicity dependence of the shape and luminos¬ 
ity of the low-mass MS. These models only go 
down to O.O83M0 and are close to the H-burning 
limit (Allard private communication), but do not 
go beyond the stellar sequence into the sub-stellar 
regime. Thus, we linearly extrapolated the BT- 
Settl models (based on the Asplund et al. 2009 so¬ 
lar abundances) for [M/H] = — 1 dex and an age of 
12 Gyr (closest to M4’s parameters) down to sub- 
stellar masses of O.O 68 M 0 , and applied distance 
and reddening as in Fig. I^for the BT-Settl model. 
The extrapolated models are plotted in Fig.[^ and 
the magnitudes around the H-burning limit are 
listed in Table As mentioned above, no sub- 
stellar models for low metallicities are currently 
available. Different models exist for a metallicity 
of [M/H] = 0.0 dex. Thus, in Fig. we compare 
our extrapolated models with these more metal- 
rich models, all for an age of 12 Gyr and scaled 
to M 4’s distance and reddening. The effect of the 
different atmospheric physics can be clearly seen 
in the shape and colour of the models. However, 
the expected end of the H-burning sequence be¬ 
tween a O.O 75 M 0 and a O.O 8 M 0 is at comparable 
magnitudes in all sub-stellar models. This gives 
us some confidence that we can use the extrapo¬ 
lated metal-poor BT-Settl models to get an esti¬ 
mate of the magnitude and colour range of the H- 
burning limit. For more exact values we will have 
to wait for the updated metal-poor models, but we 
remind the reader that the main purpose of this 
project is to provide the metal-poor benchmark 
sources and thereby hll the observational plane 
with data points that are needed to constrain the¬ 
oretical models. 

Unlike stars, BDs cannot retain their luminosi¬ 
ties via nuclear fusion. As a consequence, they 
cool with time, and a BD will be at fainter lumi- 


nosities in an old GC compared to a young BD of 
the same mass and metallicity. To get an estimate 
of this effect, we used NextGen models (Baraffe 
et al. 1997, Baraffe et al. 1998; see Fig. [^. In 
the left diagram, we plot mass against T^ff for 
various ages. The right diagram shows 1 Gyr and 
10 Gyr NextGen models for different metallicities. 
Since substellar isochrones at M 4’s low metallicity 
of [M/H]=-l dex and ages of > 10 Gyr currently 
do not exist, we also plot isochrones for a metal¬ 
licity of [M/H]=-0.5 and 0.0 dex which go down 
to lower masses. As can be seen, metallicity has a 
considerable impact on the shape of the isochrones 
as well as on the cooling time scale and hence fad¬ 
ing of low-mass, substellar objects. The [M/H]=0 
dex isochrones continue to extend with time to 
fainter magnitudes and redder colours, suggesting 
that a 0.075 Mq source fades by « 0.7 mag in 
FllOVF from 1 to 10 Gyr, and a 0.08 Mq source 
becomes fainter by « 0.2 mag. The [M/H]=-0.5 
dex isochrones suggest that a source of 0.079 Mq 
becomes fainter by « 0.7 mag, and a 0.08 Mq 
source becomes fainter by « 0.5 mag, but the 
sources also become bluer rather than redder. The 
[M/H]=-l dex 10 Gyr isochrones terminates at 
0.083 Mq. According to the models, such a metal- 
poor, low-mass source already fades by 0.3 mag in 
F\1{)W from 1 to 10 Gyr. The models suggest 
that the blue-turn is more pronounced for lower- 
metallicities, and also metal-poor sources become 
fainter compared to metal-richer sources at the 
same mass, i.e. they cool faster. 

Our optical-NlR GMDs presented in Figures 
and show that the observed MS peters out 
around F7751F « 26 mag, suggesting that we are 
approaching the end of the H-burning sequence. 
Our NIR GMD in Fig. [^suggests that the MS pe¬ 
ters out just below FWQW ~ 24 mag where the 
density of stars decreases. At fainter magnitudes, 
the WD sequence crosses the MS and the star den¬ 
sity increases again. Based on the GMDs, we es¬ 
timate our detection limits around FlTbW « 28 
mag, FllGW « 26.5 mag and Fl&GW « 27 mag. 
The detection limits are also indicated in Fig. 
with a dotted line. Gomparing the theoretical 
models in Fig. with our observed GMDs and 
taking the detection limits and the predicted H- 
burning limit into account, we conclude that we 
have reached beyond the H-burning limit in our 
NIR CMD and are probably just above or around 


this limit in our optical-NIR GMDs. 


3.4. BD candidates 

As with the best-photometry faint NIR sources, 
we visually inspected all cluster WD candidates, 
i.e., sources whose proper motion (or rather dis¬ 
placements) suggest that they are cluster mem¬ 
bers and whose position in the optical-NIR GMDs 
suggest that they are WDs. Our first photomet¬ 
ric run on the optical data did not return an opti¬ 
cal counterpart for 59 of the best-photometry faint 
NIR sources. We then over-plotted the positions of 
those 59 faint NIR sources (i.e., without an initial 
optical counterpart) on the optical F775W master 
image and inspected each position by eye. Visual 
inspection of the F777>W master image showed a 
faint optical source at or close to the location of 
the NIR source in most cases. Thus, for the second 


optical photometric run (see Sect. 2.31, the param¬ 


eter settings were relaxed so that all local maxima 
were retained, resulting in 1.5 million detections. 
Nearly all of those are just spurious detections 
(i.e., background noise or spikes in PSF wings that 
are not real stellar sources), however, a further 47 
optical counterparts to the faint best-photometry 
NIR sources were detected and thus were added 
to the initial list of optical-NIR matches. 

For the remaining twelve faint NIR sources, still 
no optical counterpart was returned. However, out 
of those, five are located on PSF streaks and two 
in the ACS WFC chip gap, so that nothing can 
be said about an optical counterpart. For four 
of the remaining five sources, visual inspection of 
the F775W master image seem to indicate an op¬ 
tical source on the position of the NIR source, but 
no photometric measurement was possible. How¬ 
ever, we used nearby stars of similar brightness 
(based on pixel counts) to estimate magnitude lim¬ 
its. One of these sources appears to have an opti¬ 
cal counterpart in the centre of our search circle, 
probably at F775W « 26 mag. This is too bright 
for a BD candidate and thus makes this source a 
WD candidate. Thus, we do not consider it fur¬ 
ther. The remaining four sources, id I to 4, with 
no optical photometry are listed in Table and 
images are shown in Fig. 13 For comparison, we 


also show images of WDs selected from the proper 
motion cleaned GMD and which have similar NIR 
magnitudes to the four NIR sources without opti¬ 
cal photometry (see Fig. 14 and Table [^. 
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Table 1: Masses and corresponding magnitudes around the H-burning limit, estimated from the extrapolated 
BT-Settl models and scaled to M4’s distance and reddening. All magnitudes are for HST/WFC3 filters and 
in the VEGAmag system. 


mass [Mq] 

FllOIU [mag] 

F160IU [mag] 

F775IU [mag] 

0.075 

25.170 

24.736 

30.105 

0.077 

24.771 

24.281 

29.327 

0.080 

24.172 

23.598 

28.160 


At the rim of our search circle for source 1, an 
optical source can be seen. This, however, is too 
far away (Ol'l away from the position of the NIR 
source) and would not agree with being a clus¬ 
ter member. In the centre of our search circle, an 
optical source just might be visible. If true, this 
source would have F775W > 28 mag, i.e., the lo¬ 
cal detection limit. For source 2, an optical coun¬ 
terpart is visible in the centre of the search circle, 
again this source would be close to the detection 
limit at F775W > 28 mag. Thus, both Source 1 
and 2 are likely (massive) BD candidates. Unfor¬ 
tunately, Source 2 is close to a saturation streak 
from a nearby bright stars in the optical image. 
The potential counterpart, however, can be clearly 
distinguished. Source 3 shows a “bright” optical 
source at the rim of our search circle, probably at 
F775W « 25.8 mag, again based on the magni¬ 
tudes of nearby stars that appear to be of similar 
brightness. However, this optical source is (/.'IS 
away from the position of the NIR source and thus 
too far away to agree with being a cluster mem¬ 
ber. A very faint optical source might just be in 
the centre of the search circle, if true, then this 
optical source would be at F775W > 28.6 mag, 
the local detection limit, making source 3 a good 
BD candidate. Just one source, source 4, does 
not show an optical counterpart at the location 
of the NIR source (i.e., within our search circle), 
and is thus our best BD candidate. Since the opti¬ 
cal photometry did not return magnitudes fainter 
than F775W = 32.9 mag, an optical counterpart 
to BD candidate 4 must be fainter than this ab¬ 
solute optical limit. Note that Source 4 appears 
somewhat extended and could probably consist of 
two or three faint sources, or possibly an extended 
object (although we expect galaxies to be much 
bluer, see Bedin et al. 2009). 

The positions of the four BD candidates are in¬ 
dicated in the CMDs in Fig. 10 Since we do not 


have an optical measurement, we provide the up¬ 
per optical magnitude limits given in Table We 
also overplot the extrapolated BT-Settl 12 Gyr 
isochrone, using the distance and reddening pa¬ 
rameters derived from the best fit in the NIR 
GMD. As can be seen, all four BD candidates are 
very close to the extension of the MS into the BD 
regime, which supports our classification of these 
sources as good BD candidates. Source 2 is blue- 
ward of the isochrone. Its position in the optical- 
NIR GMDs agrees with this source being a faint 
WD at the very bottom of the (optical) WD cool¬ 
ing sequence in M4, but it also agrees with this 
source being a VLMS star or a massive BD. In 
the NIR GMD this source is not at the very bot¬ 
tom of the WD sequence. We suggest that this 
source is a good BD candidate. Sources 1, 3 and 
4 all are very close to the MS, and their position 
in the optical-NIR GMDs does not suggest that 
these sources are WDs, but rather BDs. WDs are 
much bluer and brighter than our BD candidates. 


Figures and again show the NIR GMDs, 
but only for those NIR sources for which an 
optical counterpart had been found (black data 
points). The middle GMDs show only NIR sources 
whose optical counterparts agree with being clus¬ 
ter members, based on the displacement vector 
point diagram. We also show best-photometry 
NIR sources without optical counterparts, plotted 
in red. As mentioned in Sect. 13.11 all 177 best- 


photometry NIR sources fainter than FllQW = 
24 mag have been inspected on the FllOIU and 
the optical F775IU master images. Out of these 
faint 177 NIR sources, 165 have an optical counter¬ 
part. Out of these 165 faint NIR/optical sources, 
48 are cluster members based on their proper 
motion (displacements), and are located on the 
WD sequence. The remaining five faint best- 
photometry NIR sources which are not on PSF 
streaks or on the chip gap, and without an opti- 
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Table 2: Best-photometry NIR sources fainter than FWQW > 24 mag, and for which the optical photometry 
did not return a counterpart. The magnitudes and errors (cols. 4 and 5) are the weighted mean and error 
of the weighted mean as returned from DOLPHOT. See text for further details. 


idsD 

RA 

[h:m:s] 

DEC 
[° : ' : "] 

FllOlU 

[mag] 

F1601V 

[mag] 

F775W 

[mag] 

comment 

1 

16:23:41.701 

-26:29:17.82 

24.27 ±0.02 

23.39 ±0.11 

^28 

BD candidate 

2 

16:23:45.443 

-26:30:05.58 

25.41 ±0.05 

24.75 ±0.50 

^28 

WD/BD candidate 

3 

16:23:44.711 

-26:29:30.90 

24.36 ±0.02 

23.60 ±0.11 

> 28.6 

BD candidate 

4 

16:23:45.371 

-26:29:37.40 

26.75 ±0.16 

26.13 ±0.18 

>32.9 

BD candidate 


Table 3: Same as Tablej^but for WDs selected from the optical-NIR CMDs and with similar NIR magnitudes 
than the BD candidates. Their optical counterpart, however, is clearly detectable. The magnitudes and errors 
(cols. 4 to 6) are the weighted mean and error of the weighted mean as returned from DOLPHOT for the NIR, 
and the median-clipped mean and standard deviation for the optical magnitudes. 


iAwd 

RA 

[h:m:s] 

DEC 

[°:':"] 

FllQW 

[mag] 

F1601V 

[mag] 

F775W 

[mag] 

1 

16:23:44.354 

-26:30:26.09 

25.70 ±0.06 

25.10 ±0.24 

27.184 ±0.44 

2 

16:23:36.848 

-26:29:33.58 

26.30 ±0.11 

25.62 ±0.40 

27.618 ±0.44 

3 

16:23:36.953 

-26:29:33.81 

25.67 ±0.06 

25.14 ±0.48 

26.607 ±0.23 

4 

16:23:47.062 

-26:30:37.80 

25.68 ±0.07 

25.04 ±0.37 

26.368 ± 0.90 


cal counterpart (and hence without a proper mo¬ 
tion estimate), all make good BD candidates to 
start with. Visual inspection of these hve sources 
and estimates on the optical magnitude limit (see 
above) suggests that one is probably a WD can¬ 
didate, and the remaining four sources, listed in 
Table are good BD candidates. 

3.5. Expected number of BDs 

How many BDs can we expect? This num¬ 
ber is highly uncertain and depends on the as¬ 
sumed BD formation scenario (see Sect.[^. Fur¬ 
thermore, given our detection limits we can only 
expect to find the most massive BDs with masses 
larger than O.O 68 M 0 (based on the extrapolated 
BT-Settl models). Richer et al. (2004) derived a 
rather flat present-day mass function for M 4. Ex¬ 
trapolating towards fainter and lower-mass stars, 
they estimate that between 15 to 50 VLMSs with 
masses between 0.085 and 0.095 Mq should be in 
their held of view (GO-8679, WFPC2 data). Using 
the theoretical models (Fig. we can count the 
number of VLMSs in our held. If we only consider 
NIR sources that have an optical counterpart and 
a proper motion that suggests that they are cluster 


members, and best-photometry NIR sources with¬ 
out an optical counterpart, we hnd 23 VLMSs in 
a mass interval between 0.08 and 0.09 Mq. As¬ 
suming that the mass function is hat and that the 
slope of the mass function does not change con¬ 
siderably across the stellar/sub-stellar border, we 
can expect a similar number of BDs with masses 
between 0.070 and 0.08 Mq. However, the number 
of BDs formed per star is probably more around 
^ (see e.g. Thies et al. 2007). In this case, we can 
expect « 5 BDs. This is of course a very rough 
estimate, but it does agree with our hnding of four 
BD candidates. 

3.6. Field contamination 

How many foreground or background sources 
can we expect in our cluster CMD? We can sim¬ 
ply count the number of sources that are well 
outside the area covered by the cluster and the 
clump of held stars around AV « 0.75 pixels and 
AV « —0.5 pixels in the vector point diagram. 
We Hnd a held star density of 165 held stars per 
Apixels^. Scaling this number to the area cov¬ 
ered by our cluster stars, i.e., a displacement of 
0.1 WFC3 pixels, we hnd that we can expect 5.2 
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field stars in our cluster CMDs. 

How many foreground stars can we then ex¬ 
pect to have NIR magnitudes and colours simi¬ 
lar to our BD candidates? Selecting only field 
stars with 24 < FllQW < 27 mag and 0 < 
FllOW — F160W < 1 mag, the field star den¬ 
sity is reduced to 50 stars/Apixels^, and scaled 
to the area covered by the cluster in the vector 
point diagram , we then find that 1.6 such sources 
can be expected among our pm selected faint clus¬ 
ter members. Thus, about half of our suggested 
BD candidates might actually be foreground or 
background sources that happen to move with the 
cluster velocity across the plane of the sky. 

4. Summary and Conclusion 

We have presented the deepest NIR HST/WFC3 
study of the GC M 4 to date. The NIR data were 
proper-motion cleaned using archival deep optical 
HST/ACS (F7751T) data. Our best-photometry 
NIR CMD reveals a narrow MS delineating down 
towards the expected end of the H-burning se¬ 
quence. 177 best-photometry NIR sources fainter 
than the H-burning limit in FllGW {FllGW > 24 
mag) could be identified in our F\10W master im¬ 
age. For 165 of these faint NIR sources, an optical 
counterpart was found, 48 of these are cluster 
members according to their proper motion. All 
of these 48 faint cluster sources are on the WD 
sequence. 

We found in total five faint NIR sources for 
which the optical photometry did not return a 
measurement (and which are not on PSF streaks 
or on the chip gap). We then visually inspected 
the positions of these faint NIR sources on the op¬ 
tical images and estimated, where possible, upper 
optical magnitude limits of potential optical coun¬ 
terparts that just might be visible. One source 
is likely another WD and rejected as a BD can¬ 
didate. Based on the upper optical magnitude 
limits, we indicate the position of the remaining 
four sources in the optical-NIR CMDs. One of the 
sources (source 4 in Table , does not show an 
optical counterpart at all, which implies that its 
optical counterpart must be fainter than the ab¬ 
solute optical detection limit of FllhW > 32.9 
mag. This source appears to be somewhat ex¬ 
tended in the NIR image, which might indicate 
multiple faint sources, i.e. multiple BDs, or possi¬ 


bly a galaxy. However, its position in the CMDs 
does agree with this source being a BD. One source 
(source 2 in Table might be another WD candi¬ 
date, but its position in the optical-NIR CMD also 
agrees with this source being a massive BD or a 
VLMS star at the bottom of the MS. The remain¬ 
ing two sources also have positions that indicate 
that these sources are massive BDs. We conclude 
that we have found four good BD candidates, but 
we caution that further studies and deeper opti¬ 
cal data are necessary to confirm their status and 
cluster membership. 
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Fig. 1.— Left: Geometric-distortion corrected FllOtF master image. North is up and East to the left. The 
field of view is 136" x 123". The image is displayed on a logarithmic scale to bring out fainter sources. 
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Fig. 2.— Same as Fig. [T] but for the Fl&QW data. 
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Fig. 3.— NIR CMDs of all sources detected by DOLPHOT. Left: FllQW vs. FWQW — FIQQW, right: FIQQW 
vs. FllQW — F160W. Note that the CMDs are not selected and thus contain a large number of spurious 
detections like noise peaks and spikes around bright stars, but also faint sources that are not stars, and field 
stars that do not belong to the cluster. 
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Fig. 4.— NIR CMD that includes only best-photometry sources. The number of spurious detections should 
be minimized, however we expect that we did not include all true faint stellar sources that had originally been 
detected. The WD and BD regions have been labeled, and low-mass stellar models have been over-plotted. 
The expected end of the H-burning sequence is marked with red dashed lines and a shaded area. The magenta 
cross denotes the location of the known field BD SDSS-J125637.13-022452.4 ( |Burgasser et al. 2009[ ), scaled 
to M4’s reddening and distance. Its position in the CMD supports that our data are indeed deep enough to 
reach well into the BD zone. See text for details. 
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Fig. 5.— Top row: Vector point diagrams for all sources with counterparts with a displacement of no more 
than 2 WFC3/IR pixels. Bottom row: Optical-NIR {F775W — FllOhF) CMDs for all sources with optical 
counterparts (left); only sources with a displacement of less than 0.1 pixels, suggesting that they are cluster 
members (middle); and for the remaining field stars (right). The cluster CMD is used to select WDs, plotted 
in light blue. 
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Fig. 6.— The same as in Fig.[^ but for FllbW — F160W. 
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F775W-F110W F775W-F160W F110W-F160W 


Fig. 7.— 12 Gyr BT-Settl models (Allard et al. 2012, 2013) in HST/WFC3 filters, scaled to M4’s distance, 
reddening and metallicity. The models are extrapolated into the sub-stellar regime to 0.068 Mq. Masses are 
indicated along the sequences. The end of the H-burning sequence is estimated between 0.075 and 0.08 Mq 
and is indicated with a light-red shaded area. Detection limits in our NIR and optical data are indicated 
with dotted lines. See text for further details. 
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Fig. 8.— Comparison of our extrapolated 12 Gyr BT-Settl model with sub-stellar models at a metallicity 
of [M/H] = 0.0 dex, and with different solar abundances (AGSS2009 refers to the Asplund et al. 2009 solar 
abundance, CIFIST2011 was based on Caffau et al. 2011, ClFIST2011bc includes additional adjustments. 
See the PHOENIX webpage at https://phoenix.ens-lyon.fr/Grids/BT-Settl for more information). 
The different input physics are reflected in the different shapes of the models, however, the expected end of 
the H-burning sequence between O.Q75Mq and a O.OSMq is around FWQW « 24 mag for all models. 
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Fig. 9.— Low-mass stellar/BD models (NextGen) based on Baraffe et al. (1997, 1998), and Hausschildt et 
al. 1999. Left: Mass vs. T^ff for 1 Gyr, 2 Gyr, 3 Gyr and 10 Gyr ([M/H]=0.0 dex). Since BDs cannot 
sustain their luminosities via nuclear processes, they cool with time and become fainter. Right: Low-mass 
isochrones that show the effect of BD cooling. The dashed line is a 1 Gyr isochrone, the solid line a 10 
Gyr isochrone. Substellar models at low metallicities currently do not exist, thus isochrones for different 
metallicities are compared: [M/H]=-l dex isochrones plotted in black, [M/H]=-0.5 dex in blue, [M/H]=0 
dex in red. The [M/H]=-l dex 10 Gyr isochrone terminates at 0.083 Mq, suggesting that such a low-mass 
object already cooled and faded by 0.3 mag in FllOW. Lower masses are available for higher metallicities, 
suggesting that a 0.079 Mq source fades by 0.7 mag at a metallicity of [M/H]=-0.5 dex, about the same as 
a 0.075 Mq source at a metallicity of [M/H]=0 dex. The models suggest that lower-mass sources cool faster 
the lower the mass, and also the lower the metallicity. 
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Fig. 10.— Optical-NIR CMDs (left and middle) and NIR CMD of all sources with an optical counterpart. 
We also show the position of the four BD candidates, marked with red arrows, assuming the upper optical 
magnitude limits given in Table We also overplot the extrapolated BT-Settl 12 Gyr isochrone. As can be 
seen, the position of all four BD candidates is very close to the isochrone, suggesting that they are indeed 
BDs. See the text for more details. 
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Fig. 11.— Top row: Vector point diagram for all sources with counterparts with a displacement of no more 
than 2 WFC3/1R pixels. Bottom row: NIR CMDS {FllGW vs. FllGW - FimW) CMDs for all sources 
with optical counterparts (left); only sources with a displacement of less than 0.1 pixels, suggesting that they 
are cluster members (middle); and for the remaining held stars (right). The WDs selected from the optical- 
NIR CMDs in Figures and are plotted in light-blue, MS stars in black. Best-photometry NIR sources 
without an optical counterpart are plotted in red. As can be seen, the NIR best-photometry contains many 
VLMSs down to the expected H-burning limit, and reveals twelve sources around and below the expected 
H-burning limit {FllOW > 24 mag). Out of these, one is a likely WD candidate, whereas four sources are 
good BD candidates. See the text for further details. 
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Fig. 12.— The same as Fig. but plotted for F1601F vs. FllOlF — F1601F. 
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Fig. 13.— Zoom on the position of the 4 faint NIR sources from Table in the stacked FllQW (left) and 
FlTaW (right) images. The field of view of each image is 1.25" x 1.25" and the radius of the circles is 0.15". 
North is up and East to the left. 
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Fig. 14.— Zoom on the 4 WDs from Table in the stacked FllOVF (left) and F775W (right) images. 
Orientation and held of view is the same as in Fig.[^ The WDs have similar faint NIR magnitudes than the 
BD candidates presented in Fig. [T^ however, their optical counterpart, although faint, is clearly detectable. 
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